Linear Sigma Model at finite baryonic density and symmetry breakings 
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The linear sigma model at finite baryonic density with a massive vector field is investigated con- 
sidering that all the bosonic fields develop non zero expected classical values, eventually associated 
to condensates and corresponding to dynamical symmetry breakings which might occur in the QCD 
phase diagram. A modified equation for the classical vector field is proposed with its respective so- 
lution. Some in medium properties of the model (mainly masses) are investigated within reasonable 
prescriptions. In particular the behavior of the in medium pion and sigma masses and a particular 
way of calculating in medium coupling to baryons is investigated. A symmetry radius for finite 
baryonic densities is proposed and calculated in different ways in terms of the other variables of the 
model and these different ways of calculating it agree quite well. However, assuming that the pion 
and sigma masses go to zero close to the restoration of chiral symmetry a too high value for the 
critical density is obtained p c ~ 4.3po- 



INTRODUCTION 



Matter at high energy densities has been continuously 
investigated to provide a deepeer understanding of strong 
interactions. Experimental (relativistic and high energy) 
heavy ion collisions (r.h.i.c. and h.e.h.i.c.) provide very 
important collection of data to construct this knowledge. 
Asymptotic freedom is a key property in this program. 
With different approaches, hadronic models based on 
quantum chromodynamics (QCD) at finite energy den- 
sity have been extensively investigated and different con- 
densates are usually expected to appear at different en- 
ergy densities. In the vacuum, the lightest strong in- 
teracting particles are known to respect, approximatedly 
at least, chiral symmetry SUl{2) x SUr(2) which is ex- 
pected to be spontaneously broken down to SU(2). This 
is expected to give rise to a scalar (quark-antiquark) con- 
densate, which might have relevant effects in experimen- 
tal conditions and in realistic calculations, it rearranges 
the theory [l|,[2,S]- As a consequence, the hadronic prop- 
erties which depend on this condensate (which can be the 
order parameter or proportionally to it Q) are expected 
to vary with energy density. In particular the behavior 
of the rho vector meson has been considered as a possible 
signature of the chiral symmetry restoration via dilepton 
emission [B} although recent experimental analysis of the 
rho vector meson spectral function shows to be incom- 
patible with energy shift [7J. 

In this work the Linear Sigma Model (LSM) [8] at fi- 
nite baryonic density, ps, is investigated with a mas- 
sive classical vector field. It is based in 

Allthe 

mesons in the model are considered to develop classical 
counterparts. There is a renewed interest in the pseu- 
doscalar condensation indeed 0, HH, E2]- It provides a 
way of enhancing CP violation at finite density [l3[ . The 
exact field equations and the stability equation are trun- 
cated such as analytical solutions are obtained by consid- 
ering particular prescriptions for the stability condition. 
The numerical solutions have a self consistency although 



the so-called " full self consistency" is only achieved in a 
level of approximation and for the interactions consid- 
ered in a model, not in the complete self consistency of 
the exact realistic quantum theoretical many body prob- 
lem which still is too difficult to obtain. The trunca- 
tions in the effective action, done in the next section, are 
based in the following considerations: (1) the effective 
potential of spin zero bosons keeps the same form of that 
at the tree level calculation (i.e., quantum fluctuations 
basically rearrange the tree level model), (2) each com- 
ponent of the system, i.e. baryons/ spin zero bosons/ 
spin one fields, have nearly independent stability condi- 
tions. Hopefully this assumption might go along with 
the observation of different slope parameters and tem- 
perature freeze-out [lj| for each of the hadrons emerging 
from relativistic heavy ion collisions - and eventually dif- 
ferent contributions for the corresponding Hydrodynam- 
ics. The corresponding (dynamical) equation for each of 
the fields is satisfied. The complete numerical investi- 
gation of the results will be presented elsewhere 0, [l~5[ . 
Several properties of in medium hadrons are investigated, 
namely scalar and pseudoscalar meson masses and cou- 
plings, and their relation to the behavior of the (chiral) 
symmetry is worked out within a particular prescription 
which provides results in agreement with the usually ex- 
pected behavior. A symmetry radius is defined for the 
investigation of the symmetry behavior and its estima- 
tion (and dependence with the baryonic density) is done 
in several ways which yield very close result. In spite of 
being a quite simple model, without several degrees of 
freedom which should be relevant at high energy densi- 
ties, results show to be consistent. Assuming that pion 
and sigma masses go to zero close to the density in which 
chiral symmetry would be restored a too high critical den- 
sity is obtained. This can either signal that the model is 
too simple for describing Physics at too high densities or 
that their masses should not be expected to be zero close 
to the phase transition. 



THE LINEAR SIGMA MODEL AT FINITE p B 

The Lagrangian density of linear sigma sodel (LSM) 
with baryons, iVj(x), sigma and pions, (a,7r), covariantly 
coupled to a vector field, V^, is given by 



C = N, 



(x) (Vy^ - g s (a + i75T.Tr)) jV<(x) 



V^.Wn) + CCT + - m l v V^+ (1) 



where the covariant derivatives are given in [8| and they 
will not be completely considered in this communica- 
tion. The other terms and parameters are standard 
The introduction of a chemical potential, with an extra 
termi5£ = —NjoPchemN, is nearly equivalent to a shift 
of the classical temporal component of the vector field 
Vo coupled to the nucleons. This field however is a dy- 
namical degree of freedom (d.o.f.) and will be treated as 
such. Since the condensates (such as a =< a >) depend 
on the density, so do most of the hadronic masses. Part 
of the baryon masses are considered to come from the 
the coupling to the scalar mesonic field and part from 
an explicit mass term for the baryons in the Lagrangian: 
M* = M± g s v- 

The spin zero fields will be treated in the framework of 
the variational Gaussian approach with a truncation [3, 
9] . With the truncation of the effective potential of sigma 
and pion it can be written keeping the same form of the 
tree level effective potential. The equations for expected 
values of the sigma and pion are found accordingly, and 
shifts in these classical parts are considered due to the 
rearrangement brought by quantum fluctuations, such as 
a — > a. 

The total energy density is written in terms of the four 
variational parameters for the fields c^tt, plus baryonic 
densities and vector field variables [9j. To investigate 
the behavior the temporal component of (classical) vec- 
tor field, the total energy density is varied with respect 
to Vq, which is not quantized, instead of using its Euler- 
Lagrange equation. In this approach, Vq is found either 
by writing a reasonable (or exact) expression for ps as a 
function of Vo, as it is given below, or Vq is treated like 
a variational parameter such that one can determine a 
parametric function pb = Pb[Vo] which satisfy the equa- 
tions of mouvement and stability. For this derivation, my 
was kept constant. This second procedure yields a sort of 
variational equation for the corresponding parametric de- 
pendence. The corresponding (variational) equation for 
a constant background field component Vo can be given 
by: 



dH 
dV 



= 



gv pb 



Vo 



dV 



m 2 v Vo = 0. (2) 



where the Euler-Lagrange equation for Vo can be recov- 
ered by neglecting the derivative term above. 



The stability condition for the ground state, with 
binding energy Eq/A = H/pb < 0, can be written as 



dH 

9pb 



H 

PB 



< 0, and 



dp 2 



> 0, where po 



PB=P0 PB=PO 

is the stability density. The expressions for the energy 
density and its derivative with respect to ps is separated 
into three parts such that each component of the hadronic 
matter satisfies the stability equation above separatedly. 
With this prescription the solutions for the variational 
equations of each of the components satisfy the respective 
stability equation. The reliability of this factorization is 
not mathematically proven although some arguments for 
being reasonable are given in the Introduction and in [9| . 
The resulting equations (prescriptions) are the following: 



(i) 
(iii) 
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(3) 



dpi 



Where Hv — gvVoPB — \™vVq ''■ The complete set of 
solutions for the equations will be investigated elsewhere 
and compared to the exact numerical solutions. The vari- 
ation of the sigma and pion masses arise from the cor- 
responding classical fields. In the vacuum p 2 can go to 
zero as long as tt 2 — > satisfying the Goldstone theorem 
when c — ► in the Lagrangian. 



Densities, coupling constants and masses 

The baryon fields, which depend on the bosonic fields 
through the Dirac equation coupled to the mesons, are 
quantized in terms of creation and annihilation operators. 
The baryonic degrees of freedom sum up into the densi- 
ties: baryonic (ps), scalar (ps) and pseudo-scalar (p ps ) 
densities. These quantitites will not be explicitely evalu- 
ated here although they are partially used below 0, [l(| • 
The energy density due to the fermions (antifermions) 
(Ef j) and the density of baryons (antibaryons) (Pb,b) 
can be written, in the leading order, in terms of (I) their 
momenta for each kind of baryons (i), up to the last oc- 
cupied level with momentum fc^, and of (2) the classical 
vector field as 0, [l6| : 



E 



(/,/) ~ V*) 



^ d 3 k ( 2 -^~ 



(M;+^ +] _ ) )+V-o(Vb-2g ( i +i _ ) ) 



B,B 



— f k * d 3 k 



2(M>+E\ + _ ) ) 

\_ (4) 

In these expressions E± — gVo ± \/k 2 + (M*) 2 are the 
eigenvalues of the corresponding Dirac equation [l(| E3] • 
These expressions still correspond to an approximation 
and show deviations from the Fermi liquid picture. 

The scalar and pseudoscalar densities, which appear in 
the equations of cr, 7f, can be expanded in terms of the 
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scalar and pseudoscalar condensates, for example, as: 



for which a 2 +ir < v 2 , at ps = Po, it follows respectively: 



PS 



PPS 



(o) 
Ps 



a (1) 
-Ps 



erf. 



&vac 

n\Pps + F \Pps 

M 



(2) 



\k \Pps 



a 3 
a 3 

vac 



(3) 

Ps 



■o(a), 



(5) 

where the coefficients pg'pg are obtained from the ex- 
pressions calculated with the solution of the correspond- 
ing Dirac equation. They are such that ps = p ps = 
when a = a vac and re = 0, with higher order terms are 
indicated by o(|tt|) and o(a). By substituting these ex- 
pressions into the variational equations for the respective 
condensates Q it is found that the terms proportional to 
a and tv can yield contributions to the in medium masses 
of sigma and pions. The terms proportional to a 3 and 
•fr can produce (effective) contributions for the coupling 

constant, i.e., A — * A* = A ± j§— ~ A ± p PS - These two 
corrections for the in medium effective coupling constant 
may also be different from each other, eventually lead- 
ing to different interactions of the pion and sigma in the 
baryonic medium. 

The second equation of prescriptions ^ can produce 
the same solution of the equation f2}. The third of 
prescriptions §3§ can define a symmetry radius in the 
medium: 



-2 



- 2 



(6) 



In this expression C is a constant to be determined from 
the parameters of the model. In the vacuum: a 2 = v 2 = 
f 2 as discussed above. Modifications in the equation |[3l- 
{Hi)) will produce different dependences on the baryonic 
density. A more general symmetry radius, corresponding 
to particular modifications of the corresponding differen- 
tial equation due to diverse couplings for example, might 
be written as: (a 2 + -k 2 — v 2 ) = (D + CpB c ) 7 where 
D, C, c, 7 are constants to be related to the parameters 
of the model. 

One way of calculating C is found by assuming that 
this symmetry radius is valid over a range of baryonic 
densities. This is a crude approximation because, heavier 
hadrons as well as quark and gluon d.o.f. are expected to 
be relevant for high energy densities. In the high density 
when chiral symmetry should be restored: a — 7r — > 0. 
This critical density is written as p c = upo,. Thus in this 
point: C = Tv 2 y-^^- At the saturation density (ps) 
the expression for C can be written as: 




U\2 




(7) 



Four values are considered: (i) u = 2, (ii) u — 3, (iii) u — 
3.5 and (iv) u = 4. Considering the branch of solutions 



Vz 2 



* 2 



* 2 



I'D 



0.54w (i), 
0.68?) (iii), 



Vt 2 - 
\fd 2 



0.65w (ii), 
~ 0.71v (iv). 
(8) 



Other solutions are not presented. Since the squared 
value 7f 2 is a scalar which appear very often in the expres- 
sions, it may be that /* ~ \/ b 2 + ■n 2 , i.e., a pion clas- 
sical field could be responsible for modifications in the 
pion decay constant and consequently measurable, even 
if competing with other effects. The topological Skyrme 
model can provide some argument in favor of such in- 
terpretation for a classical pion field inside hadrons 17 1. 
These expressions may be therefore useful for relating 
descriptions of different ranges of the matter phase dia- 
gram. 

The values obtained for C from estimates (jHJ) are re- 
spectively given by: 

(7~±0.41fnr5 (i), (7~±0.33fm~5 
C~ ±0.30fm"5 (iii), C~±0.28fm"^ (iv). 



(9) 



Another way of estimating C is shown by considering 
the meson masses in the medium. With the expressions 
for meson masses in terms of the classical fields and v 
the symmetry radius can be written as: C^fps = 

^ ) ) 2 + (^ ) ) 3 at 



MPT? - K ac ) 2 ), where (^ } ) 2 



a given density pg. In these expressions the coupling 
A was also kept constant (and positive) and c = 0, in 
the Lagrangian term. Two possible behaviors are ob- 
tained in this picture for the restoration of chiral sym- 
metry: the sum of these masses may decrease or increase 
depending on the sign of C^fps. For po = 0.15fm~ 3 
and C ~ — 0.15fm - ^ the above expression yields ap- 
proximated values (pt) 2 (po) — (1 ± 0.53) Pt(pb = 0). 
If one considers that the pion and sigma masses disap- 
pear close to the chiral symmetry restoration point (i.e., 
if (p*) 2 — » 0), with the values above we obtain that 
p c ~ 4.3/90. Seemingly it is a too high baryonic density 
and the reasons are quite apparent. Firstly, as empha- 
sized above, the present work only takes into account the 
light sector of hadrons and it does not consider quark and 
gluon degrees of freedom. Furthermore, it is a controver- 
sial subject whether pion and sigma masses (two point 
Green's functions) should be expected to be so close to 
zero (as it was assumed to obtain such high value for 
the critical density) close to (and at the) deconfincment 
critical point. 

These ways of calculating C provide crude (but in- 
teresting and curious) estimations. The corresponding 
in medium hadron properties are qualitatively in agree- 
ment with other estimations Q ■ The inclusion of other 
relevant d.o.f. will be presented elsewhere as well as a 
corresponding calculation at finite temperature. 
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Summary and Conclusions 



Munhoz and G. Krein. 



In this work some aspects of the Linear Sigma Model 
at finite baryonic density were investigated with a mas- 
sive classical vector field, based in [tj |T(|. All the mesons 
in the model were considered to develop classical coun- 
terparts. In part this is due to independent new inves- 
tigations on pseudoscalar condensates which have shown 
a renewed interest in the pseudoscalar condensation in- 
deed 0, EH, [IH . The exact field equations and the stabil- 
ity equation were truncated for obtaining analytical so- 
lutions which capture the expected behavior of the sys- 
tem. These solutions have a self consistency although 
the so-called " full self consistency" is only achieved in a 
level of approximation and for the interactions consid- 
ered in a model, not in the complete self consistency of 
the exact realistic quantum theoretical many body prob- 
lem which still is too difficult to obtain. The trunca- 
tions in the effective action, done in the next section, are 
based in the following considerations: (1) the effective 
potential of spin zero bosons keeps the same form of that 
at the tree level calculation (i.e., quantum fluctuations 
basically rearrange the tree level model), (2) each com- 
ponent of the system, i.e. baryons/ spin zero bosons/ 
spin one fields, have nearly independent stability con- 
ditions. Hopefully this assumption might go along with 
the observation of different slope parameters and temper- 
ature freeze-out for each of the hadrons emerging from 
relativistic heavy ion collisions - and eventually different 
contributions for the corresponding hydrodynamic. The 
corresponding (dynamical) equation for each of the fields 
are satisfied. The complete numerical investigation of the 
results will be presented elsewhere 0, HH . Several prop- 
erties of in medium hadrons were investigated, namely 
scalar and pseudoscalar meson masses and couplings, and 
their relation to the behavior of the (chiral) symmetry is 
worked out within a particular prescription which pro- 
vides results in agreement with the expected behavior. A 
symmetry radius was defined for the investigation of the 
symmetry properties and its estimation (and dependence 
with the baryonic density) is done in several ways with 
fair agreement, in spite of being a quite simple model, 
without several degrees of freedom which should be rele- 
vant at high energy densities. Related aspects to matter- 
antimatter asymmetry in relativistic heavy ion collisions 
and in the Early Universe will be discussed and investi- 
gated elsewhere 

Hi- 
Acknowledgements 

This work was partially supported by FAPESP, Brazil. 
F.L.B. thanks brief discussions with L.McLerran, R. 
Rapp, D. Zschiesche, H. Stoecker, R. Pisarski, M. 



[1] S. Weinberg, The Quantum Theory of Fields -II, Cam- 
bridge University Press, (1996). 

[2] F.L. Braghin, Doctoral thesis, (1996) unpublished. 

[3] P.M. Stevenson, B. Alles and R. Tarrach, Phys. Rev. 
D 35 2407(1987). F.L. Braghin, Phys. Rev.D 57, 3548 
(1998); Phys. Rev. D 64, 125001 (2001). F.L. Braghin 
and F.S. Navarra, Phys. Lett. B 508, (2001) 243. 

[4] H. Umezawa, Advanced Field Theory, AIP Press, (1995). 
J. Stern, Nucl. Phys B-Proc. Suppl. 64, 232 (1 998). 

[5] T. Hatsuda and T. Kunihfro, |nucl-th/0112027] G. Chan- 
fray, M. Ericson, P.A.M. Guichon, Phys. Rev. C 63, 
055202 (2001). D. Roder,J. Ruppert, D.H. Rischke, Phys. 
Rev. D 68 016003 (2003). J. Wambach, Prog, in Part, 
and Nuc. Phys. 50 615 (2003). R. Rapp and J. Wambach, 
Adv. Nucl. Phys. 25, 1 (2000). D. Jido, T. Hatsuda, T. 
Kunihiro, Phys. Rev. D63, 011091 (2001). J. Delorme, 
G. Chanfray, M. Ericson, Nuc. Phys. A 603, 239 (1996). 
L.S. Celenza, H.S. Wang, CM. Shakin, Phys. Rev. C 63, 
025209 (2001). GE. Brown and M. Rho, Phys. Rep. 269, 
333 (1996). N. Kaiser, S. Fritsch, W. Weise, Nuc. Phys. 
A 697, 255 (2002). P.K. Panda and F.L. Braghin, Phys. 
Rev. C 66, 055207 (2002). 

[6] K. Su zuki et al, Ph ys. Rev. Lett. 92, 072302 
(2004); |nucl-ex/0211023| K. Yokokawa, T Hatsuda, A. 
Hayashigaki, T. Kunihiro, |hep-ph/020"4163"j 

[7] Damjanovic for the NA60 Collaboration, arXiv: 
|nucl-ex/0 609026 (2006). 

[8] S. Gasiorowicz and D.A. Geffen Rev. of Mod. Phys. 41, 
531 (1969). 

[9] F.L. Braghin, pre-print IF-1570/2001; |hep-ph/0206107l 

IF-2006. F.L. Braghin, in HADRON Physics 2002, ed. 

by C. Vasconcellos et al World Scientific (2002). 
[10] F.L. Braghin, Int. Journ. of Mod. Phys. D 13, 1267 

(2004). F.L. Braghin, accepted to be published in Int. 

Journ. of Mod. Phys. D. 
[11] H. Shen, S. Tamenaga, H. Toki, |nucl-th/0409057| 
[12] A.B. Migdal et al, Phys. Rep. 192 179 (1990). 
[13] See, for example, Y. Mori, K. Saito Phys. Lett. B 552, 

21 (2003). 

[14] For example in: F. Grassi, Braz. J. Phys. 35 (2005) 52 
and references therein. HADRON- RANP Joint Meeting, 
RJ, Brazil, 2004, ed. by M. Bracco et al, AIP Conference 
Proceedings 739, APS (2004). 

[15] F.L. Braghin, pre-print IF-USP (2005). 

[16] F.L. Braghin, works presented in: ISMD 2006, 
Paraty, RJ, Brazil, September 2006; Physics in Col- 
lisions, Buzios, RJ, Brazil, July 2006 -e-Proceedings: 
|arXiv:h ep-ph/0609242 18 th International Conference 
on Few Body, Santos, SP Brazil, August 2006. 

[17] For example: T.H.R. Skyrme, Proc. Roy. Soc. Ser A 
260, 127 (1961). GS. Adkins, in Chiral Solitons ed. by 
GS. Adkins, Add. Wesley, 1988. F.L. Braghin, LP Cav- 
alcante, Phys. Rev. C 67 065207 (2003) and references 
therein. 



4 



